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Abstract 

Aim of this paper is the qualitative analysis of the solution of a 
fH ■ boundary value problem for a third-order non linear parabolic equation 

which describes several dissipative models. When the source term is 
linear, the problem is explictly solved by means of a Fourier series with 
properties of rapid convergence. In the non linear case, appropriate 
estimates of this series allow to deduce the asymptotic behaviour of 
the solution. 
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1 Introduction 

We refer to the non linear equation 

(1.1) CsU = F{x,t,U,Ux,Ut) 

where Ce is the third - order parabolic operator: 

(1.2) C, = d,,{edt + c^)-dt{dt + a). 

In (1.1)-(1.2) e,a,c are positive constants and F = F{x,t,u,Ux,ut) is a 
prefixed function. 
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The equation (1.1) characterizes the evolution of several dissipative mod- 
els such as the motions of viscoelastic fluids or solids, the sound propagation 
in viscous gases, the heat condution at low temperature and some Joscphson 
effects in superconductivity when dissipative causes are not neglegible (for 
details see, e.g, [6] [7]). 

The analysis of several initial - boundary problems related to the equa- 
tion (1.1) has been discussed in many papers (see, for example [10] and its 
references) . 

In [6], the fundamental solution K(x,t) of the operator has been 
determined explictly by means of convolutions of Bessel function, and so, 
when F = f{x,t) is linear, the initial value problem has been discussed. 
Further, the linear strip problem has been solved in [8] by means of a Green 
function G connected with K{x,t). However, it seems difficult to deduce an 
exhaustive asymptotic analysis by means of this solution. 

In this paper, both linear and non linear cases are analyzed. The Green 
function of the linear strip problem is determined by Fourier series which 
is characterized by properties of rapid convergence and allows to establish 
an exponential decrease of the solution. Moreover, the non - linear problem 
is reduced to an integral equation with kernel G and, by means of suitable 
properties of G, the asymptotic analysis of the solutions is achieved. 



2 Some typical problems in superconductivity 

As an example of physical models related to the equation (1.1), we will 
consider a typical case of non linear phenomenon concerning the Josephson 
effects in Superconductivity. 

Let if = (p{x, t) be the phase difference of the wave functions related to 
the two superconductors, and let 7 be the normalized current bias; when 
F = sirup — ^, the equation (1.1) gives the perturbed Sine- Gordon equation 
[1]: 

(2.1) eipxxt + ^xx - ^tt - aipt = sin (/? - 7. 

The terms ecpxxt and aipt characterize the dissipative normal electron 
current flow along and across the junction, respectively, and they represent 
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the perturbations with respect to the classic Sine Gordon equation [11, 14, 
16]. When the surface resistence is neglegiblc [4, 12, 13], then e is vanishing 
and a singular perturbation problem for the equation (1.1) could be appear. 
As for the cocfRcicnt a of (2.1), it depends on the shunt conductance [5] 
and generally one has a « 1 [3, 13, 14]. However, if the resistence of the 
junction is so low to short completely the capacitance, the case a » 1 
arises [1, 17, 18]. 

If I is the normalized length of junction, a first boundary problem refers 
to the phase - boundary specifications for x = and x = I [5]. Moreover, 
according to the junction geometry, other typical boundary value problems 
could be considered. For instance, in an overlap geometry [1, 14], if r] is the 
normalized measure of the external magnetic field applied in the plane of 
the junction [4], it has: 



Instead, if an 'open-circuit' is examined, then boundary conditions (2.2) 

assume the homogeneous form [3], while, for the electrodynamics of a Joseph- 
son junction in two spatial dimensions, other conditions on (0,1) can be 
considered [15]. 

As an example, we will analyze the problem related to the phase bound- 
ary specifications but the analysis can be applied to the other problems. 



3 The linear strip problem 

Let Z,T be arbitrary positive constants and let 



(2.2) 



iPx{^,t) = ifx{l,t) = T]. 



9, = {{x,t):Q<x<l, < i < T}. 



In the linear case, the strip problem is: 



(3.1) 



dxxi^ut + c^u) - dtiut + au) = f{x, t), {x, t) G O, 
u{x,0) = go{x), ut{x,0) = gi{x), xe[0,l], 
u{0, t) = 0, u{l, t) = 0, 0<t<T. 
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As it's well known, the problem (3.1) can be solved explicitly by means 
of the Fourier method, as soon as the Green function for the operator (3.1)i 
has been determined. For this let : 

TtTT 1 / 

(3.2) 7„ = — , bn = cjn, K = -^{a + e-fl), ujn = \J - b^; 
and 

(3.3) Hn{t) = —e-'^"'^sinh{cOnt)- 

By standard techniques, the Green function related to (3.1) can be given 
the form: 

2 °° 

(3.4) G{x,t,^) = j^Hnit) sin'jnC sin-fnX. 

n=l 

whose properties will be discussed in the next section 4. 
Then, if we consider the functions: 

(3.5) ug^= f 9i{0G{x,^,t)d^ i = 0,l 

Jo 

and 

(3.6) uf = G*f=['dT f /(^, r) G(x, ^, i - r) d^, 

JQ Jo 

the formal solution of (3.1) is given by: 

(3.7) u{x, t) = Ug^ + {dt + a- edxx)ug^ - Uf{x, t). 
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4 Properties of the G function 



At first we will analyze the convergence and the asymptotic properties 
of the series (3.4) which defines G. We will prove that the order of the terms 
of G is at least and the function G is exponentially vanishing as t ^ oo. 
For this we put: 

(4.1) P= , + a{l/.Y ^ " = 2 ' ^ = 

and we will prove the following lemma: 

Lemma 4.1-Whatever the constants a,e,c^ may be in 3?+, the function 
G{x, ^, t) defined in (3.4) and all its time derivatives are continuous functions 
in Q. Moreover, everywhere in J7, it results: 



f0 G 

(4.2) \G{x,(,,t)\ < Me-^\ \^\ < Nj e'^^ j G N 



where M,Nj are constants depending on a, e, c^. 

Proof. When = ae, the operator can be reduced to the wave 
operator. If < ae, for all n > 1, one has /i„ > bn- So, if A; is an arbitrary 
constant such that c^/ae < k < 1, the Cauchy inequality assures that: 

(4.3) ^ <\{sjl + a) = hn. 
Then, if X„ = (bn/hn)^ , one has: 

(4.4) (1 _ A;)i < ^ = (1 - Xn)-^ < 1 - XJ2. 
Further, for all ri > 1, it results: 
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(4.5) ^^>W''' =(«-«/2)"' 
and hence, by (4.4) and (4.5) i, it follows: 

(4.6) (,-t{hn-OJn) < g-pt_ 

Moreover, the estimates (4.4), (4.5)2, (4.6) imply that the terms iJ„ 
defined by (3.3) verify: 



(1-A;)-V2 e-P* 
g — a/2 



(4.7) \Hn\ < \ 2-, Vn > 1 



and hence, (4.2)i follows. As for (4.2)2, we observe: 



6^ 2c2 

(4.8) /in - u;n = , " < — , Vn > 1 

hn + iOn £ 

and by means of standard computations one obtains (4.2)2. 
Consider now the case > ae, and let 



(4.9) iVi,2 = -[lT (l-a£/c2)V2]. 

If /c is a positive constant less than one, let iVi_25 be the lowest integers 
such that 



(4.10) N*< Ni, N* >N2; , Nk > + {1 - aek/c'^)^/% 

£7rv k 

It suffices to analyse only the hyperbolic terms with n > Simi- 
larly to (4.4) for all n > Nk, it results : 



(4.11) COn >hn{l- kf'^. 
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So (4.5)2, (4.6), (4.8) and (4.11) imply the estimates (4.1) also when (? > as. 

Obviously, the lemma holds also when ATi 2 are integers. In this case, 
the constants M and Nj could depend on t. □ 

As it is known, as for the x-differentiation of Fourier's series like (3.4), 
attention must be paid to convergence problems. For this, we will consider 
x-derivatives of the operator [edt + c^)G instead of G and Gt- 

Lemma 4.2.- Whatever a, e, may be , the junction G{x, ^, t) defined in 
(3.4) is such that: 

(4.12) \eGt + c^G\ < e-^\ 

(4.13) \d,,{eGt+^G)\ <A^ e-^\ 

where Aq and Ai are constants depending on a,e,(P. 

Proof. As for the hyperbolic terms in we consider Xn = (^)^ < 1 
and put: 

(4.14) ^^ = hn{-l + ^/T^^). 
It results: 

(4.15) eHn + c^Hn = 77^e-^"*{(c2 + ev?n)e""* - [c^ - s{hn + u;n)]e-""*}, 

ZUn 

where Taylor's formula assures that: 

(4.16) c +Sipn = — — 2 - ^hnXn - — / 7572 dy. 

a + s% 16 Jo (1 - yr''' 

Further, it is possible to prove that, for all n, it is: 
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(4.17) < 2cl/s7r) 
while, for every h> and n > a{l + h), one has: 

(4.18) l""" dy<lx'^ {{h + lf/[h{h + 2)]V2 _ 1}. 

So, estimates (4.17) and (4.18) allow to obtain: 

(4.19) \c^ + eipn\ <\ (5aV4 + ai/8), 

where ai = a'^Hh + lf/{h^ + 2hf/^ - 1]. 

The estimate of Lemma 4.1 together with (4.19) show that the terms of the 
serie related to the operator (4.12) have order at least of n~^. So, it can be 
differentiated term by terms with respect to x and the estimate (4.13) can 
be deduced. □ 

As solution of the equation CeV = we will mean a continuous function 
v{x,t) which has continuous the derivatives vt,vu, dxx{£Vt + (?v) and these 
derivatives verify the equation. 
So, we are able to prove the following theorem: 

Theorem 4.1- The function G{x,t) defined in (3.4)?s a solution of the 
equation 

(4.20) C,G = dUeGt + (?G) - dt{Gt + aG) = 0. 

Proof. The uniformly convergence proved in Lemma 4.1 and Lemma 
4.2 allow to deduce that: 



2 °° 

{dt+a)Gt = J ^{[-bl+{2hn-a)hn]Hn+{a-2hn)e''''''^ coshuntjsmjn^smjnX, 

n=l 

(4.21) 



8 



and 



dxx{£dt+c^)G = y 51 [i^llhn-bl)Hn-e-ile coshwnt] sin7„^sin7„a;, 

n=l 

(4.22) 

and hence, (4.20) foUows.D 

5 Properties of the convolution 

To obtain the exphcit solution of the hnear strip problem, the properties 
of the convolution of function G with the data must be analyzed. 
For this, let g{x) be a continuous function on (0, 1) and let consider: 

(5.1) ug{x,t) = tgiO G{x,^,t) d^. 

Jo 

(5.2) u*g{x,t) ={dt + a-edxx)ug{x,t). 

So the following lemma holds: 

Lemma 5.1- If g{x) is a C^{0,1), then the function Ug defined by (5.1) 
is a solution (n.4) of the equation = ^ such that: 

(5.3) limug{x,t) = lim dtUg{x, t) = g{x), 

i ^ i ^ 

uniformly for all x G (0, 1). 

Proof. Lemma (4.1)-(4.2) and continuity of g assure that function (5.1) 
and the partial derivatives required by the solution converge absolutely for 
all {x,t) G O. Hence, by theorem 4.1, one deduces: 

(5.4) dxx{edt + c^)ug = f g{0 dU^dt + ^)G{x, ^, t) d^ = dt{dt + a)ug. 

Jo 
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Besides, the hypoteses on the function g and (4.2), imply (5.3)i. 
Moreover, one has: 



(5.5) Gt = — > Hnit) siwynX, 

TT of n 

^ n=l 



and hence: 



2 H it) 

(5.6) dtUg = ^ — — [5(0cos7„^]o sin7„a; + 



TT ' — n 

n=l 



TT Jo n 

n=l 

So, if r]{x) is the Heaviside function, from (5.6), one obtains: 



(5.7) nni9ttx, = f[5(Z)-5(0)]+5(0)- /^'(0[^(e-^) + f-l]^^e = 5(a^)-° 

t-^O I Jq I 



Lemma 5.2- Ifg{x) G C^{0,1) with g^'\0) = g'^'\l) = 0(i = 1,2,3), then 
the function u* defined in (5.2) is a solution (n.4) of the equation = 
such that: 



(5.8) \imu*g{x,t) = g Ihn dtu*g{x,t) = 

uniformly for all x G (0, /). 

Proof. The hypoteses on g{x) allow to have 

(5.9) OxxUg{x,t) = g"{^) G{x,^,t) d^ = Ug»{x,t). 
So, by Lemma 5.1 and expression (5.9), one obtains : 
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(5.10) dt{dt + a)u*g = dt{dt + a)[{dt + a)ug - eu = dxx{sdt + c ) 



Moreover, it results: 



(5.11) 



dtu* = dfUg + adtUg - edtu » = c^u » , 



and hence, (5.3) i implies (5.8)2- 
Finally, owing to (5.3)-(5.9), one has: 



(5.12) 



lim u* = lim[{dt + a)ug — £u n] = g{x).U 



6 Explicit solution of the strip linear problem 



Let us consider the homogeneous case . ^Prom Lemma 5.1 and Lemma 
5.2 the following result is obtained: 

Theorem 6.1- // the initial data gi{x), and go{x) verify the hypotheses 
of Lemma 5.1 and Lemma 5.2, then the function: 



(6.1) u{x, t) = Ug^ + {dt + a- edxx)ugo 

represents a solution (n.4) of the homogeneous strip problem: 



3.2) 



dxx{eut + c^u) - dt{ut + au) = 0, {x, t) G Q, 

u{x,0) = go{x), ut{x,0) = gi{x), x£[0,l], 
. u(0,t) = 0, u{l,t) = 0, te{0,T]. 



As for the non-homogeneous problem, let: 
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(6.3) Uf{x,t) = f'dT f fi^,T) G{x,^,t-T) d^, 

Jo Jo 

the following theorem holds: 

Theorem 6.2- // the function f{x, t) is a continuous function in Q with 
continuous derivative with respect to x, then the function Uf represents a 
solution (n.4) of the non-homogeneous strip problem: 

dxx{sdt + (?)u - dt{dt + a)u = -f{x, t), {x, t) G 

(6.4) i u{x,0) = 0, dtu{x,0) = 0, xe[0,l], 

, u{0,t) = 0, u{l,t) = 0, te{o,T]. 

Proof. By (3.4) one has: 



(6.5) 

and hence: 



lini 'u/(a;, t) = 0, 



(6.6) dtUfix,t) = f dT f f{i,T) Gt{x,i,t-T) d^. 

Jo Jo 

Moreover, as in (5.7), it is possible to prove that: 



(6.7) lim /'/(C,T) Gtix,^,t-T) d^ = f{x,t), 

and one obtains: 



(6.8) d^Uf =f{x,t)+ f dr f f{i,T) Gu{x,i,t-T) d^. 

Jo Jo 

So that. Lemma 4.2, Theorem 4.1 and (6.4)-(6.6) assure that: 
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d^^{edt+c^)uf = j dr j fi^,T)d^^{eGt+c^G)d^ = dtidt+a)uf-f{x,t) 



Jo Jo 



(6.9) 



Finally, owing to (6.3)-(6.6) and Lemma 4.1, if Bi (i=l,2) are two positive 
constants, one has: 



from which (6.4)2 follow. □ 

Moreover, the uniqueness is a consequence of the energy-method (see, 
for example, [9]) and we have: 

Theorem 6.3- When the source term f{x,t) satisfies theorem 6.2, and 
the initial data {go,gi) satisfy theorem 6.1, then the function 

(6.11) n(x, t) = Ug^ + {dt + a- edxx)uga - uj 

is the unique solution (n.4) of the linear non-homogeneous strip problem 



At first we will consider the linear problem (3.1). Obviously, the asymp- 
totic properties of the solution depend on the behaviour of the source term 
/ = f{x,t). For instance, one has: 

Theorem 7.1- When the source term f{x,t) satisfies the condition: 



(6.10) 



Uf\ < Bi{l - e-^^)- \dtuf\ < ^2(1 - e-^^), 



(3.i)n 



7 Asymptotic properties 



(7.1) 



f{x,t)\ < h 



1 



{h, k,a = const > 0) 



{k + ty+' 



-a 



then the solution of (3.1) is vanishing as t — >■ 00, at least as t 
Moreover, when 
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(7.2) \fix,t)\ < C e-^* (C, 5 = const > 0), 

one has: 



(7.3) \u{x,t)\ < ke ^'^ S* = min{l3,S}, k = const. 

Proof. Referring to the explicit solution (6.11) of the problem (3.1), we 
observe that, owing to the Lemma 4.1, the terms Ug^ and Ug^ defined by 
(5.1) arc exponentially vanishing as t — t- oo . Moreover, the convolution Uf 
given by (3.6) can be estimated by (4.2) too and the hypotheses on f{x, t) 
imply the statement. □ 

Consider now the non linear problem: 



(7.4) 



suxxt + c^Uxx - utt - aut = F{x, t, u), 
u{x,0) = go, ut{x,(}) = gi, 
u{0,t)=0, u{l,t) = 0. 



When one applies the results of section 6 (see 6.3), one obtains the 
following integral equation: 



(7.5) u{x, t) = f gi {OG{x, ^, t)d^ + (9* + a - ed^,) f go{OG{x, ^, t)d^+ 

Jo Jo 

f dr f Gix, i,t-T) F{i, T, u{i, T))di, 
Jo Jo 

to which the previous estimates can be applied. 

For instance, as for the perturbed Sine Gordon equation (2.1), one has: 

(7.6) \F{x,t,u)\ = \sinu + ^\ < 71, 71 > 

and so the Lemma 4.1 implies that the solution of the related strip problem 
is bounded for all t. 
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Another example of applications could be given by exponentially vanish- 
ing non linear source F(x,t,u(x,t)) such that (see [2]): 



(7.7) \F{x,t,u{x,t))\ < const, e'"* (/^ > 0) 

As consequence of (7.5) and (4.2), the solution of (7.4) will vanish expo- 
nentially, too. 
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